Variational calculation of the ppK system based on chiral SU(3) dynamics 
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The ppK~ system, as a prototype for possible quasibound K nuclei, is investigated using a 
variational approach. Several versions of energy dependent effective KN interactions derived from 
chiral SU(3) dynamics are employed as input, together with a realistic NN potential (Avl8). Taking 
into account theoretical uncertainties in the extrapolations below the KN threshold, we find that 
the antikaonic dibaryon ppK~ is not deeply bound. With the driving s-wave KN interaction the 
resulting total binding energy is B(ppK~) = 20 ± 3 MeV and the mesonic decay width involving 
KN — > ttY is expected to be in the range 40 - 70 MeV. Properties of this quasibound ppK~ system 
(such as density distributions of nucleons and antikaon) are discussed. The A(1405), as an I — 
quasi-bound state of K and a nucleon, appears to survive in the ppK~ cluster. Estimates are 
given for the influence of p-wave KN interactions and for the width from two-nucleon absorption 
(KNN — * YN) processes. With inclusion of these effects and dispersive corrections from absorption, 
the ppK~ binding energy is expected to be in the range 20 - 40 MeV, while the total decay width 
can reach 100 MeV but with large theoretical uncertainties. 

PACS numbers: Valid PACS appear here 



I. INTRODUCTION 

In the context of low-energy QCD with Nf = 3 quark 
flavors, the study of possible antikaon- nuclear quasi- 
bound states is a topic of great current interest. Sponta- 
neously broken chiral SU(3) x 5/7(3) symmetry, together 
with explicit symmetry breaking by the non-zero quark 
masses, basically determines the leading couplings be- 
tween the low-mass pseudoscalar meson octet (Nambu- 
Goldstone bosons in the chiral limit) and the octet of 
the ground state baryons. In particular, the Tomozawa- 
Weinberg chiral low-energy theorem implies that the 
driving KN interaction in the isospin 1 = channel is 
strongly attractive. Likewise, the 1 = 7rS interaction 
is attractive. The coupling between these KN and 7rE 
channels is the prime feature governing the subthreshold 
extrapolation of the KN interaction. A detailed knowl- 
edge of this subthreshold interaction is required when 
exploring the possible existence of bound A"-nuclear clus- 
ters. 

The quest for strong binding and dropping antikaon 
masses in a nuclear medium has a long history which 
originated in early discussions of kaon condensation in 
dense matter [1, 2] and continued over the years [3-6] 
with increasing levels of refinement [7]. The topic has 
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recently been revived when the existence of long-lived, 
deeply bound ^-nuclear states was suggested using a 
simple potential model [8]. It was argued that if the 
-ftT-nuclear binding is sufficiently strong to generate such 
systems below 7rE threshold, their width could indeed be 
small [8, 9]. Experiments performed in search for such 
states [10-12] have so far not given conclusive answers. 
While the search continues with more detailed analy- 
ses [13-17], "non-exotic" final state interaction scenarios 
cannot be ruled out as interpretation of the data [18]. 

An important prototype system for these considera- 
tions is ppK~ , the simplest antikaon-nuclear cluster. It 
has recently been investigated using three-body (Fad- 
deev) methods [19, 20] and variational approaches [21- 
24] . Reaction studies [25] have also been performed deal- 
ing with the actual formation of ppK~. The Faddeev 
and variational calculations predict a total ppK~ bind- 
ing energy in a range B ~ 50 - 70 MeV, together with an 
estimate of the KNN -> nYN decay width, r ~ 50 - 100 
MeV, depending on details of the interactions used. 

The key issue in any such calculation is the (model de- 
pendent) extrapolation of the KN interaction into the 
region far below threshold. Its predictive power is so 
far limited by the persistent lack of accurate constraints 
from data. It is in this perspective that we perform the 
present variational calculation using a subthreshold effec- 
tive KN interaction systematically derived from chiral 
SU(3) coupled-channel dynamics. 

Apart from the constraints provided by KN thresh- 
old data and low-energy cross sections, the only piece 
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of information about the interaction below KN thresh- 
old is the 7rS mass spectrum which is dominated by the 
A(1405) resonance. The subthreshold extrapolation of 
the effective KN interaction has recently been investi- 
gated in detail [26] from the viewpoint of chiral SU(3) 
dynamics [27-31]. In this approach the off-shell KN am- 
plitude below threshold is governed by the strong attrac- 
tion in the KN and 7rE channels, and by the dynamical 
coupling of these channels. These couplings are again de- 
termined by the Tomozawa- Weinberg chiral low-energy 
theorem. Their structure shares features with the early 
pioneering coupled- channel model [32] that used vector 
meson exchange interactions (see also Ref. [33]). 

Most chiral SU(3) based calculations agree that the 
A(1405) resonance structure, with its nominal position at 
1405 MeV as seen in the 7rS mass spectrum, is actually 
shifted to about 1420 MeV in the KN amplitude as a con- 
sequence of coupled-channel dynamics. This implies that 
the effective single-channel KN interaction is substan- 
tially weaker than anticipated in the simple phenomcno- 
logical potential used previously in Rcfs. [8, 21]. In those 
phcnomenological studies, the local, energy-independent 
potential was adjusted interpreting the A(1405) directly 
as a KN bound state, identifying its binding energy by 
the location of the maximum observed in the 7rE spec- 
trum, but ignoring strong coupled-channel effects. 

In this paper we perform a variational ppK~ calcula- 
tion employing the new effective KN potential derived 
from chiral coupled-channel dynamics [26] , together with 
a realistic NN potential. This calculation is supposed to 
be complementary to the Faddeev approach with chi- 
ral SU(3) constraints [20]. The variational calculation 
gives easy access to the wave function of the bound state 
so that valuable information about the structure of the 
ppK~ cluster can be extracted, whereas the elimina- 
tion of the 7rS channel is required and the width of the 
state can only be estimated perturbatively. The Faddeev 
calculation has, in turn, the advantage that the decay 
width of the quasibound state is computed consistently 
in the coupled-channel framework. Both methods there- 
fore have their virtues and limitations which need to be 
discussed in comparison. 

The present work extends and improves our previous 
studies [24] in several directions, including further re- 
finements in the NN interaction, computation of density 
distributions, an evaluation of effects from p-wave KN 
interactions and an estimate of the KNN — ► YN ab- 
sorptive width. The paper is structured as follows. The 
variational framework and formalism are developed in 
Sec. II. The derivation of the effective KN interaction 
based on the chiral SU(3) coupled-channel approach is 
briefly summarized and discussed in Sec. III. Results of 
our ppK~ calculations are presented in Sec. IV followed 
by a summary and conclusions in Sec. V. 



II. FORMALISM 

A. Model wave function 

In search for the energetically most favorable KNN 
configuration, the present variational investigation fo- 
cuses on the ppK~ system with spin and parity J 71 " = 0~ 
and isospin (T,T Z ) = (1/2,1/2), where the parity as- 
signment includes the intrinsic parity of the antikaon. 
We prepare the following two-component variational trial 
wave function: 



\*)=tf- 1 [\* + )+C\*-)], (1) 

where A/" -1 is a normalization factor and C is a mixing 
coefficient. The components |$+) and |3>_) have the form 



1*+) 
|*_) 



<S> + { ri ,r 2 ,r R ) \S N = 0) 

x [[NN] Tn=1 K] t=1/ ^ Tz=1/2 

$-(n,r 2 ,r R ) \S N = 0) 

[I nn ]t„=oK] t=1/2iTm=1/2 



(2) 
(3) 



with spatial wave functions $±{r\,r 2 ,r R ) multiplied by 
the spin state vector of the two nucleons and the isospin 
state vector of the total KNN system, respectively. In 
both components, the spin of the NN pair is assumed 
to be zero (Sn = 0). The large component |<f> + ) has 
isospin (T/v = 1) of the two nucleons corresponding to the 
dominant ppK~ configuration, with inclusion of pnK° 
through charge exchange. The admixture of the compo- 
nent [<!>_} (Tjv = 0), with zero isospin in the NN sector, 
can occur through a combination of I = and 1=1 
KN interactions and turns out to be small, typically less 
than 5%. Both components have the same total isospin 
(T=l/2). 

The detailed ansatz for the spatial wave functions is 
chosen as follows: 

$±(ri,r 2 ,r R ) = F N (n) F N (r 2 ) F K (r R )G(r 1: r 2 ) 
x [Hi{r u r R )H 2 {r 2 ,r R ) 

±H 2 (r 1 ,r R )H 1 (r 2 ,r R )} .(4) 

Here Fjv(rj) {i = 1,2) and F R (r R ) are trial functions 
describing the localization of the nucleons and the kaon, 
respectively. Their forms are assumed to be single Gaus- 
sians: 



F N {ri) = cxp[-^r 2 2 ] , F K {r R ) = exp[- 7 r|] . (5) 

The /j, and 7 in Eq. (5) are treated as independent pa- 
rameters. Correlations between the two nucleons, and 
between the K and each nucleon, are described by cor- 
responding correlation functions: 



N N 



G(r 1 ,r 2 ) = l-^g„exp -A„ (n - r 2 ) 

71=1 

N K 

H a (ri,r R ) = 1 + ^ /i a>rl exp -v n {n-r R ) 



71=1 



(6) 



(7) 
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The NN correlation function, G(r\,r2), is prepared to 
account for the strong short-distance repulsion in the NN 
interaction which keeps the two nucleons apart. The KN 
correlation functions, H a (ri,r R ), are given the flexibil- 
ity to adjust themselves appropriately to the attractive 
antikaon-nucleon interaction. 

One notes that the spatial wave functions 
<&±(ri, r 2 , r R ), Eq. (4), are even or odd under ex- 
change of the two nucleons: 

<&±{r 2 ,ri,r R ) = ± $±(n, r 2 , r R ) . (8) 

The NN pair in 1$+) (1$-)) is thus in a singlet-even 
(singlet-odd) state. The spatial wave function is rotation- 
ally symmetric, i.e. the total orbital angular momentum 
is L = 0. 

The trial wave function has the following real-valued 
variational parameters: C in Eq.(l), \x and 7 in Eq. (5), 
{flVi,A„} (n = 1,...,N N ) in Eq. (6), and {h a . n ,v n } (n = 
1, ...,Nk) in Eq. (7). The range parameters of the Gaus- 
sians in the NN and KN correlation functions are orga- 
nized as 



The energy of the center-of-mass motion, 



X n — X 



1 • 



Ai 



JVjv-1 



v\ ■ 



VNk 



N K -1 



(9) 



independently for each of the NN or KN sectors. Var- 
ious combinations of parameter sets for the Gaussian 
ranges, /1, 7, (Ai, Xn n ) and (fi, vn k ) have been tried. 
For each combination of the range parameters, we find 
a set of parameters C and {g n ; h a . n } which minimizes 
the expectation value of the total Hamiltonian, using the 
Simplex method [34]. This determines the ppK~ bound 
state of minimal energy, if existent. 

A remark should be added concerning the treatment 
of the center-of-mass (CM) motion. Given the indepen- 
dence of the variational parameters [i and 7 in Eq. (5), 
the CM wave function cannot simply be separated by fac- 
torization in the present model. A complete separation is 
possible only in the special case j/jj, = ttik/Mn, where 
Mjy (mx) is the nucleon (kaon) mass. In the actual cal- 
culations we have confirmed that this relation between 
7 and n turns out to be satisfied quite accurately even 
though it has not been imposed from the beginning. The 
variational procedure favors indeed a wave function in 
which the CM motion factorizes as it should. One typ- 
ically finds [i ~ 0.2 fm~ 2 and 7 ~ 0.1 fm~ 2 , using the 
KN interaction based on chiral SU(3) dynamics and de- 
scribed in Sec. III. 



B. Hamiltonian 

The Hamiltonian used in the present study is of the 
form 

H = f + V NN + RcV RN -f CM • (10) 
Here T is the total kinetic energy: 

f = fcp| + iL. (ii) 



Tcm = 



(P1+P2+PR) 
2 (2M N + m K ) 



(12) 



is subtracted. 

As a realistic nucleon-nucleon interaction Vjviv we 
choose the Argonne vl8 potential (Avl8) [35]. Since the 
total spin of the two nucleons is restricted to zero as ex- 
plained in the previous section, the tensor, LS and (LS) 2 
potentials do not contribute. We thus employ the central, 
L 2 and spin-spin parts of the Avl8 potential: 

Vnn = ]T P(X) 

X= 1 E, 1 

x [v c x (r) + v x 2 (r)L 2 + v s x s (r) a, ■ «r 2 ] , (13) 

where P(X) is a projection operator onto the singlet- 
even ( 1 i?) or singlet-odd ( 1 0) state, and L is the orbital 
angular momentum operator for the relative coordinate 
between two nucleons. The central and L 2 potentials 
[v x (r) and v x 2 (r)] are identified with the phenomeno- 
logical short- and intermediate-range parts of the Avl8 
potential, Eq. (20) in Ref. [35]. The long-range spin-spin 
term comes from one-pion-exchange, Eq. (17) in Ref. [35]. 
We ignore the electromagnetic part of Avl8. Both 1 E 
and O potentials are taken into account since our model 
wave function includes both types of NN states. In the 
present study the r-dependence of each of the potential 
terms is well fitted by a series of Gaussians (see Ap- 
pendix A). The most pronounced feature is the strong 
short-distance repulsion in the singlct-cvcn central po- 
tential. 

The energy-dependent effective s-wave KN interaction 
Vk N is represented as 



Vrn = HKN 1 )+v(KN 2 ) , 



v(KN) = Pi(KN) 

■ kn(Vs) cxp [-{r RN /a s f] 



1=0,1 



x v{> 



(14) 



(15) 



where r RN 



2M N 2m K 



'k ~ r N for each of the two nucleons, A^i 
and N2, and Pi(KN) is the isospin projection operator 
for the KN pair. In the present work, the radial de- 
pendence of the KN potential is assumed to be a single 
Gaussian form with range parameter a s . The interaction 
strength v KN (y/s) is a function of the (off-shell) center- 
of-mass energy variable y/s of the KN system. This in- 
teraction, extrapolated into the subthreshold region, is 
a key issue in the present paper and will be specified in 
greater detail in a separate Sec. III. 



C. Calculational procedure 

The energy dependence of the KN interaction requires 
a self-consistent variational procedure to minimize the 
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energy of the KNN system. This is done in the same 
way as in our previous work [22] . 

We introduce an auxiliary (non-observable) antikaon 
"binding energy" Bk to control the CM energy y/s of 
the KN subsystem within the ppK~ cluster. This Bk is 
defined as 



-B K = (tf - (*|#iv|*> 



(16) 



where Hn is the nucleonic part of the Hamiltonian, 

Hn = Tn + Vnn — Tcm,n, (17) 

rf, PI+PI rf (P1+P2) 2 , . 

I N — —^T7 ) I CM,N = J77 • (i-O) 



2M N 



The relation between the KN two-body energy y/s and 
Bk within the three-body system is not a priori fixed. 
In general, 



y/s = M N + m K - 77 B K 



(19) 



where 77 is a parameter describing the balance of the an- 
tikaon energy between the two nucleons of the KNN 
three-body system. One expects 1/2 < r\ < 1. The 
upper limit (7/ = 1) corresponds to the case in which 
the antikaon field collectively surrounds the two nucleons, 
a situation encountered in the limit of static (infinitely 
heavy) nucleon sources. In the lower limit (77 = 1/2) the 
antikaon energy is split symmetrically half-and-half be- 
tween the two nucleons. We investigate both cases and 
label them "Type I" and "Type II" , respectively: 

Type I : yfs = M N + m K -B K , (20) 
Type II : y/s = M N + m K - B K /2 . (21) 

The calculations then proceed as follows. First, as- 
sume Bk to be some trial starting value, Bjp. Given a 
relation between Bk and y/s, the strength of the KN po- 
tential is now fixed, and the Hamiltonian is determined. 
Then the variational calculation is performed to find the 
state of minimal energy. Given that state, a new an- 
tikaon binding energy, B^\ is evaluated with the wave 
function so obtained. Then one examines whether B^ 
coincides with Bjp. If not, a different starting value B 



K 



is chosen and the procedure is repeated until B^) = B^ 



'K 



is satisfied at an acceptable level of accuracy. 

The KN potential is in general complex. In order to 
perform the variational calculation of the energy and the 
"bound state" the real part of the potential, Re V K n> 
is used as a starting point. The decay width T of that 
state is then calculated perturbatively by taking the ex- 
pectation value of the imaginary part of the KN poten- 
tial: 



r 



-2 (y\ImV KN \y 



(22) 



At this stage the width T represents the mcsonic two- 
body decay channels (KN — > 7rE, it A) within the KNN 
three-body system. The non-mesonic absorption width 
for KNN -> EN, AN will be treated separately. 



III. EFFECTIVE KN INTERACTION 

A. KN potential based on chiral SU(3) dynamics 

Here we discuss the effective KN potential developed 
in Rcf. [26]. This potential has been systematically con- 
structed using chiral SU(3) coupled- channel calculations 
which successfully describe S = — 1 meson-baryon scat- 
tering and the properties of the dynamically generated 
A(1405) resonance. The formulation of the coupled- 
channel approach is briefly sketched in Appendix B. 
Starting from this coupled-channel framework, a complex 
and energy-dependent interaction kernel V eS is derived in 
the single KN channel, such that the full coupled-channel 
s-wave KN scattering amplitude is exactly reproduced 
by solving the single-channel equation with V oS . Then 
this effective interaction kernel is approximated by an 
equivalent local KN potential used in the Schrodinger 
equation (with reduced KN mass /i), 



d 2 u(r) 



+ 2fi [U(r, y/s) + B] u(r) = . (23) 



The potential U(r, y/s) is expressed by the effective inter- 
action kernel V eS (y/s) together with a normalized spatial 
distribution g(r), 



U{r,yfs) = 



g(r) M N 

2U) y/s 



V eS (yTs), 



(24) 



where r is the relative coordinate of the KN system and 
<D = [s 2 — (Mjf — m 2 K ) 2 ]/ (4s 3//2 ) is the reduced energy of 
the KN two-body system. 

An important observation in Ref. [26] and earlier work 
is that the 7rS diagonal coupling is strong enough to gen- 
erate a resonance in the single (elastic) 7rS channel. In 
the strongly interacting KN <-» 7rE system, the A(1405) 
appears as a KN bound state embedded in the resonant 
7rS continuum. The experimentally observed broad spec- 
trum in the 7rE channel, with its maximum at 1405 MeV, 
is dominated by the inherent 7rS interaction, while the 
resonant structure in the KN amplitude, governed by 
the inherent KN interaction, actually appears around 
1420 MeV. The KN "binding energy" commonly asso- 
ciated with the A(1405) is therefore not 27 MeV but 
only less than half of this naive estimate. This small 
KN binding energy implies a weaker attractive poten- 
tial, only about half as strong as the phenomenological 
potential of Refs. [8, 21]. One should nonetheless note 
that both phenomenological and chiral potentials repro- 
duce the existing experimental data around threshold. 
Their qualitatively different subthreshold extrapolations 
result from the fact that the coupled-channels framework, 
constrained by the chiral effective Lagrangian, induces 
off-shell dynamics which is very different from the purely 
phenomenological approach. 

Here we choose the spatial distribution of the potential 
as a Gaussian form 



U(r,y/s) = w^ w (Vs) cxp [-(r/a s f 



(25) 
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where a s is the range parameter. The strength of the 
approximate local potential is related to the interaction 
kernel V cS {y/s) defined in Ref. [26] as 



KN 



M N 



2tt 3 / 2 al u^fs 



(26) 



The potential so obtained is complex and energy depen- 
dent, reflecting the elimination of the other (mainly irY>) 
channels. The center-of-mass energy y/s is related to the 
binding energy B of the KN two-body system B as 



B 



fs + M N + m K ■ 



We choose the range parameter a s such that the reso- 
nance structure is reproduced at the position predicted 
by the full chiral dynamics calculation. The scattering 
amplitudes for both 1 = and 1=1 around KN thresh- 
old are also well reproduced by this potential. 

However, the effective single- channel KN effective in- 
teraction V is generally non-local. Naive translation of 
V eS into an approximate local form does not guarantee 
that this local potential (we refer to it as "uncorrected" ) 
reproduces the KN amplitude of the full coupled-channel 
calculation over a wide range of subthreshold energies. It 
is indeed found that simple extrapolation of the local po- 
tential (26) to the deep subthreshold region, y/s < 1400 
MeV, significantly overestimates the scattering ampli- 
tude in comparison with that of the original coupled- 
channel approach [26] . Compensation of this deficiency 
requires modifying the strength of the real part of the 
potential and introducing extra energy dependence. The 
strengths of these "corrected" potentials in the 1 = and 
1=1 channels are parametrized by polynomials as 



KN 



= Kr + K Itl s 1 ' 2 + Kr 2 S + Kr 3 s 



8 3/2 



(27) 



The coefficients Ki i are given in Ref. [26]. The strengths 
of the "corrected" and "uncorrected" potentials at r = 
are shown in Fig. 1 , based on the chiral model of Ref. [37] . 
For the subsequent variational three-body calculation we 
always use the "corrected" potentials. 

In order to estimate theoretical uncertainties we use 
altogether four different variants of chiral dynamics 
calculations from which we derive equivalent energy- 
dependent local potentials, ORB [36], HNJH [37], 
BNW [38], and BMN [39]. All these models reproduce 
the total cross sections for elastic and inelastic K~p scat- 
tering, threshold branching ratios, and the 7rS mass spec- 
trum associated with the A(1405). Differences among 
those models mainly stem from the lack of accurate data 
for the 7rS spectrum and from differences in details of 
the fitting procedures. The values of the range param- 
eters for these models are shown in the second row in 
Table I. Additional uncertainties, concerning higher or- 
der terms in the interaction kernel derived from chiral 
SU(3) dynamics, are estimated to be about 20 % around 
y/s ~ 1360 MeV, based on the systematic study of such 
higher order corrections in Ref. [38]. 





-200 



> 



> 



-500 



■r---- 1 . i 1 i 
_ (a) 1=0 

- 


1 1 

- 


uncorrected 



i 


1320 1360 1400 
s 1/2 [MeV] 


1440 


.__ 1 1 1 1 1 


1 1 
1=1 


' — corrected 




uncorrected 





1320 



1360 



1400 



1440 



s 1/2 [MeV] 



FIG. 1: (Color online) Strength of the "corrected" (solid lines) 
and "uncorrected" (dotted lines) potentials at r = with 
HNJH model [37]. The real parts are shown as solid/dotted 
lines. Imaginary parts are given as dashed lines, (a): 1 = 
channel, (b): 1 = 1 channel. 



B. Structure of the A(1405) 

Let us now examine the two-body KN system with 
1 = where the A(1405) is generated dynamically below 
the KN threshold. In order to study the structure of 
the A(1405), we first solve the Schrodinger equation (23) 
with the real part of the potential, Re[f/(r, \/s)]. This 
treatment appears to be justified by the relatively small 
imaginary part of the potential as seen in Fig. 1. Cor- 
rections from the dispersive shift of the binding energy 
induced by the imaginary part of the potential will be 
estimated as we move along. 

The self-consistent results for the binding energies are 
summarized in Table I together with the root mean dis- 
tance (y/ (r 2 )) between antikaon and nucleon which form 
the A(1405). These results are produced with the "cor- 
rected" potentials, but the "uncorrected" ones give es- 
sentially the same output since the difference in strength 
between the "corrected" and "uncorrected" potentials is 
small in the energy region relevant to the KN bound 
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TABLE I: Range parameters a a of the local KN potential, 
the self-consistent KN binding energy B, and the root mean 
distance \J (r 2 ) between antikaon and nucleon. The energies 
B± and B2 are determined, respectively, by the zero of the 
real part and the maximum of the imaginary part of the full 
amplitude with complex potential (see text). 





ORB [36] 


HNJH [37] 


BNW [38] 


BMN [39] 


a a [fm] 


0.52 


0.47 


0.51 


0.41 


B [MeV] 


11.8 


11.5 


9.97 


13.3 


Vtr*) [fm] 


1.87 


1.86 


1.99 


1.72 


Bi [MeV] 


16.5 


15.5 


17.0 


16.7 


B 2 [MeV] 


17.8 


16.9 


19.8 


18.9 



state (y/s ~ 1420 MeV). The typical KN binding energy 
and root mean distance found with the chiral potential, 

B ~ 12 MeV , VT^>" ~ 1-9 fm , (28) 

should be compared with the results of the p henomeno- 
logical model [21]: B = 27 MeV and VT^I = 1-36 fm. 
The small binding obtained with the present potential 
is related to the strong 7r£ interaction in chiral dynam- 
ics, as discussed in Ref. [26]. Since the binding energy is 
smaller, the size of the bound state becomes correspond- 
ingly larger than that of the phenomenological model. 

It is instructive to recall the study of electromagnetic 
properties of the A(1405) in chiral dynamics [40], where a 
relatively large electric mean squared radius has been re- 
ported. Although the electric mean squared radius is not 
directly comparable to the mean distance of antikaon and 
nucleon, qualitative agreement of these independent size 
estimates gives support to the chiral effective potential 
introduced here, while the radius in Ref. [40] was com- 
puted with the full chiral coupled-channel amplitudes. 

C. Dispersive effects induced by the imaginary 
part of the potential 

The variational calculation of the ppK~ three-body 
problem has the disadvantage (unlike the Faddeev ap- 
proach) that only the real part of the complex KN poten- 
tial can be handled as input whereas the imaginary part 
must be treated perturbatively. It is therefore manda- 
tory to estimate the systematic uncertainties caused by 
this limitation. 

In the KN two-body case, the scattering amplitude 
can easily be obtained by solving Eq. (23) with the full 
complex potential. Effects of the dispersive shift on the 
binding energy, induced by the imaginary part of this 
potential, can then be examined by comparing the full 
result (with complex potential) to the one obtained us- 
ing only the real part of the potential, see Eq.(28). The 



energy, B\, of the subthreshold KN state generated by 
the complex potential can be deduced from the zero of 
the real part of the KN scattering amplitude. For com- 
parison we also check the position of the maximum, B 2 , 
of the imaginary part of that amplitude, which coincides 
with B\ if the non-resonant background is small. 

The calculated values of B\ (and B 2 ) are listed, for 
all models considered, in Table I. Comparing B\ with 
the binding energies B found with the real parts of the 
potentials, we conclude that the dispersive effects may 
increase the binding energy of the KN two-body quasi- 
bound state by 

AB ~ 3 - 7 MeV, (29) 

or AB ~ 6 - 10 MeV when comparing B 2 and B. We 
keep this shift AB in mind for a discussion of systematic 
uncertainties when we now turn to the KNN three-body 
system. 



IV. RESULTS 

This section presents the results of our variational cal- 
culations of the quasibound ppK~ system. Four vari- 
ants of effective KN interactions have been employed: 
ORB, HNJH, BNW and BMN as explained in the previ- 
ous section. These interactions are translated into local, 
energy-dependent potentials. In all cases we have used 
the "corrected" versions of the potentials which prop- 
erly reproduce the KN scattering amplitude computed 
with full chiral coupled-channel dynamics. In the self- 
consistent treatment of the energy dependence of the in- 
teraction, both options, "Type I" [Eq. (20)] and "Type 
II" [Eq. (21)] for the relationship between y/s and Bk 
have been tried. With these different choices for on-shell 
equivalent potentials and different options for handling 
off-shell energies within the three-body system, we are in 
a position to give rough estimates of the uncertainties as- 
sociated with the required subthreshold extrapolations. 
In all calculations, the convergence of the Gaussian ex- 
pansions (6) and (7) has been checked and found satis- 
factory with Nn — Nk = 5. 



A. Total binding energy and decay width 

The results of self-consistent solutions for total ppK~ 
binding energies are summarized together with the 
KNN -> ttYN decay widths in Table II and depicted 
in Fig. 2. In all present calculations the binding ener- 
gies turn out to be quite modest. Deeply bound, nar- 
row ppK~ states are not seen to develop. The varia- 
tional calculations, with KeVg; N as previously specified 
and lmVx N treated according to Eq.(22), predicts total 
binding energies and decay widths in the range 

B(ppK~) - 20 ± 3 MeV and T - 40 - 70 MeV 
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TABLE II: Results of variational calculations for ppK~ bind- 
ing energies and (KNN — > ttYN) decay widths F based 
on four versions of chiral effective KN potentials. "ORB" 
etc. indicate the chiral SU(3) model used in each case. The 
"corrected" equivalent local potentials have been adopted 
throughout, with ranges a s given in Table I. "Type I (II)" 
refers to the ansatz for the relationship between y^s and Bk 
(see text). 
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BNW 


BMN 


"Type I" 












B(ppK-) 


[MeV] 


18.9 


16.9 


18.1 


16.6 


r 


[MeV] 


53.1 


47.0 


60.4 


38.8 


"Type II" 
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[MeV] 
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r 


[MeV] 
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53.1 
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FIG. 2: Distribution of total binding energy and mesonic 
decay width. The models "ORB", "HNJH", "BNW" and 
"BMN" are shown with symbols square, circle, triangle and 
diamond, respectively. Closed (Open) symbols indicate the 
"Type I" ("Type II") ansatz. 



for all cases studied. The "Type II" ansatz favors slightly 
stronger binding (20 - 23 MeV) than the "Type I" op- 
tion (16 - 19 MeV). Different potentials produce binding 
energy variations within only about two MeV for each 
given "Type I" or "Type II" set. 

It is instructive to examine the detailed decomposi- 
tion of the total ppK~ energy into kinetic and poten- 
tial energies of the antikaon and two-nucleon subcom- 
ponents. These sets of numbers are given in Tables III 
and IV for the four variants of chiral SU(3) based models 
mentioned previously. Shown (separately for the "Type 
I" and "Type II" options) are the total kinetic energy 
Ekin — (^\T— Tcm\^), the nuclear part of the kinetic en- 
ergy, T nuc = (*|Tjv-Tcm,jv|*), the antikaon binding en- 
ergy Bk defined in Eq. (16), the expectation value of the 
if-nuclear potential energy, V(KN) = (<f|Re V KN \^), 



and the contribution from the nucleon-nucleon interac- 
tion, V(NN) = (*|t>ivjv|*}- 

Several interesting observations can be made. First, 
the kaon "binding energy" Bk is in the range 40 - 50 
MeV for all cases studied. Note, however, once again that 
Bk is not an observable in the KNN three-body system. 
In fact the observable total ppK~ binding energy is less 
than half of Bk- The nucleons arc the slow movers in 
the quasibound compound: their kinetic energies add up 
to only about 40 MeV, while the total kinetic energy in 
the three-body cluster is typically more than three times 
larger and thus carried predominantly by the antikaon 
floating between the two slowly moving, heavy nucleons. 
The antikaon's potential energy roughly cancels the total 
kinetic energy, leaving room for the nucleon-nucleon in- 
teraction to bind the system which, as an isolated proton- 
proton pair, would be unbound. 

The admixture of the isospin-zero (Tjv = 0) compo- 
nent |<f>_) [Eq. (3)] is typically about 4%. It originates 
from the coupling matrix element (&+\Vkn\®-) which 
is proportional to the difference of / = and 1=1 KN 
interactions. Although small, this admixture helps bind- 
ing the ppK~ system: without the \&-) component, the 
total binding energy would decrease by about 5-7 MeV 
as seen from the lines denoted "HNJH^" in Tables III 
and IV. 

The spin-spin and L 2 parts of the NN potential have 
an influence on the \$~) admixture. The contribution of 
these terms to V(NN) is small and attractive. It tends to 
reduce the mixing of the |$_) component into the total 
wave function. This is seen in the last lines of Tables III 
and IV denoted by "HNJH*" where the spin-spin and L 2 
potentials are turned off and only the central part of the 
NN potential is active. The spin-spin and L 2 terms were 
not considered in our previous Ref. [24]. 



B. Structure of the ppK cluster 

Consider now the more detailed structure and charac- 
teristic sizes of the ppK~ system as found in the present 
calculations. For this purpose we introduce density dis- 
tributions of the NN and KN pairs in the quasibound 
compound as functions of the respective nucleon-nucleon 
and antikaon-nucleon distances, as follows: 

Pnn{x) = (*| S 3 ( ri -r 2 -sc)|tf) , (30) 

p KN (x) = \ £ SHrK-n-x)^) , (31) 

i=l,2 

where r\ t 2 refer to the two nucleons and tk to the an- 
tikaon. Both densities are understood as being normal- 
ized to one. The projected density distributions for KN 
pairs with specific isospin I = 0, 1 are 

Pkn( x ) = 
(*| \ Pi(KN t )6 3 (r R -r t -x)\*) 



(32) 



= 1,2 
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TABLE III: Detailed compilation of ppK~ results calculated for chiral SU(3) models ORB, HNJH, BNW and BMN (see 
text) with the "Type I" ansatz for the relation between antikaon binding energy Bk and y/s in the off-shell KN two-body 
subsystem. The listing includes: the total binding energy B(ppK~) and decay width Y for KNN — > nYN; the total kinetic 
energy, Ekin = (^\T — Tcm I*}, the nuclear part of the kinetic energy, T nuc — (*|Tjv — Tcm,n\^)\ the antikaon binding energy 
Bk as defined in Eq.(16); the if-nuclear potential energy, V(K N) — (^Re Vjj ]v | , I'}; the contribution from the nucleon-nucleon 
interaction, V(NN) — (^l Vjvjvj'I'), with all entries given in MeV. Lower part of table: percentage P(Tn = 0) of the 1$-} 
admixture to the total wave function; r.m.s. distances Tijvjv and Rkn between the two nucleons, and between antikaon and a 
nucleon; r.m.s. distances between antikaon and nucleon, R KN (I = 0, 1), for isospins I = 0, 1 of the KN pair (all given in fm). 
Lines denoted HNJFL are calculated switching off the admixture of the |$-} component. Lines HNJH* are results obtained 
without the spin-spin and L 2 potentials of the NN potential. (Both these cases refer to the KN potential derived from the 
HNJH model). 





B(ppK~) r 






T 

J- nuc 


B K 






V(KN) 


V(NN) 


ORB 


18.9 53.1 




125.9 


38.3 


40.2 






-127.7 


-17.0 


HNJH 


16.9 47.0 




129.5 


38.1 


38.9 






-130.4 


-16.2 


BNW 


18.1 60.4 




124.9 


37.4 


39.2 






-126.7 


-16.3 


BMN 


16.6 38.8 




141.5 


39.9 


40.8 






-142.5 


-15.7 


HNJHt 


12.0 44.8 




115.1 


37.0 


31.0 






-109.0 


-18.2 


HNJH* 


15.9 47.1 




129.6 


37.3 


38.9 






-131.2 


-14.3 






P(T N = 0) 


Rnn 




Rkn 


Rkn {1 — 


0) 




Rkn (I = 1) 


ORB 


3.4 % 


2.15 




1.93 




1.79 






2.28 


HNJH 


3.8 % 


2.21 




1.97 




1.82 






2.33 


BNW 


4.1 % 


2.20 




1.97 




1.81 






2.35 


BMN 


3.8 % 


2.23 




2.00 




1.86 






2.34 


HNJH* 


% 


2.13 




2.01 




2.01 






2.01 


HNJH* 


4.5 % 


2.26 




2.00 




1.83 






2.39 


TABLE IV: 


Same as Table III, but using the 


"Type II" ansatz for the relation between antikaon binding energy Bk and 


in the off-shell KN two-body subsystem. 




















B{ppK-) r 






Tnuc 


B K 






V(KN) 


V(NN) 


ORB 


22.7 64.2 




136.0 


40.9 


46.1 






-141.2 


-17.6 


HNJH 


20.8 58.3 




141.0 


40.9 


45.1 






-145.2 


-16.6 


BNW 


18.1 71.4 




132.0 


39.3 


43.2 






-136.0 


-16.7 


BMN 


21.7 53.1 




158.4 


43.8 


49.3 






-163.9 


-16.2 


HNJHt 


13.8 51.8 




121.3 


38.7 


33.9 






-116.5 


-18.5 


HNJH* 


19.8 58.6 




141.5 


40.3 


45.2 






-146.4 


-14.9 






P(T N = 0) 


Rnn 




Rkn 


Rkn (I — 


0) 




Rkn (I = 1) 


ORB 


3.7 % 


2.10 




1.85 




1.71 






2.21 


HNJH 


4.1 % 


2.15 




1.89 




1.73 






2.26 


BNW 


4.1 % 


2.16 




1.91 




1.75 






2.29 


BMN 


4.1 % 


2.15 




1.89 




1.73 






2.25 


HNJHt 


% 


2.09 




1.96 




1.96 






1.96 


HNJH* 


4.8 % 


2.19 




1.91 




1.74 






2.31 



with Ni = J d 3 x p^jylx). These average distances are 
summarized in the lower parts of Tables III and IV. The 
relatively weak binding of the system implies rather large 
NN distances, typically around Rnn — 2.2 fm, while the 
average KN distances are slightly smaller, R^ N — 1-9—2 
fm. When looked at separately in the 1 = and 1=1 
channels of the KN subsystem, the significantly stronger 
attraction in the 1 = component drags the antikaon 
closer to the nucleon than in the 1=1 component. 

The different density distributions, Eqs. (30-32), are 



with the isospin projectors Pi(KNi). 

With these distributions one can determine mean- 
square distances 

R 2 nn = J d 3 x x 2 pnn(x) , (33) 

R\ N = J d 3 xx 2 PRN (x) , (34) 

Rr n (I) = NJ 1 J d'xx^^^x) , (35) 
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shown in Figs. 3 and 4. These densities depend only on 
the absolute value of the relative coordinate x. The two- 
body densities are plotted as functions of r = \x\ and 
normalized as 4n J drr 2 p(r) = 1. The two-nucleon dis- 
tribution pmn shows the pronounced effect of the short- 
distance repulsive core of the NN interaction. The maxi- 
mum NN density reached in this weakly bound system is 
about 0.03 fm -3 , a small fraction of typical bulk nuclear 
densities and only twice as large as the maximum proton- 
neutron density in the even more dilute, very weakly 
bound deuteron. 

The KN density distribution has its maximum at zero 
distance between the antikaon and each nucleon. It re- 
flects the strong KN attraction in the 1 = channel, 
whereas the 1=1 KN density is small. The right panel 
in Fig. 4 shows a comparison between densities of the 
1 = KN pair in the ppK~ cluster and the one form- 
ing the A(1405) as a two-body KN quasibound state. 
Both densities are properly normalized for the compar- 
ison. The conclusion to be drawn from this picture is 
that the A(1405) stays essentially intact in the ppK~ sys- 
tem which appears to behave much like a weakly bound, 
short-lived pA(1405) compound. Again, the 1=1 com- 
ponent of the total KN density is small. A detailed anal- 
ysis of the expectation values of the squared angular mo- 
mentum, L 2 ^ N in the ppK~ system shows that it is close 
to zero in / = and close to two in I = 1 KN configura- 
tions, indicating as expected the dominant s-wave in the 
1 = channel. The survival of the A(1405) in the three- 
body cluster is qualitatively consistent with the result of 
the phenomenological potential [21]. 



C. Comparison with Faddeev results 

The small ppK~ binding energy of only about 20 MeV 
found in the present variational approach appears to be 
inconsistent with results of a Faddeev calculation using a 
chiral interaction [20] which predicts about 80 MeV bind- 
ing. Here we discuss possible reasons for this difference. 

An advantage of the Faddeev treatment is its capabil- 
ity to treat the three-body dynamics in the 7rEA^ chan- 
nel. The variational approach, on the other hand, works 
with an effective KN interaction after eliminating the 
7rS channel. While the attractive irY, two-body inter- 
actions and their coupled-channel effects arc nonetheless 
accounted for as part of the complex and energy depen- 
dent KN potential, additional attraction may indeed be 
generated by the ttT,N three-body dynamics treated ex- 
plicitly in the Faddeev approach. 

Secondly, there are significant differences in the sub- 
threshold extrapolations of the two-body KN interac- 
tion. When solving Faddeev equations, a one-term sep- 
arable form has been used to approximate the KN in- 
teraction [20]. While this interaction agrees with ours 
in its on-shell properties around KN threshold, it has 
been pointed out in Ref. [26] that using such a sepa- 
rable approximation gives stronger attraction at lower 



(subthreshold) energies, as an artifact of the regular- 
ization procedure. This is also evident by comparison 
of subthreshold extrapolations with early works on chi- 
ral SU(3) coupled-channel dynamics which had adopted 
separable forms for the KN and 7r£ interactions [27]. 
In contrast, the KN potential in the present investiga- 
tion is constructed using dimensional regularization of 
loops in the Bethe-Salpeter equation 1 (see Appendix B 
for a further detailed assessment of different regulariza- 
tion schemes). The resulting subthreshold amplitudes 
have significantly smaller real parts in comparison with 
those using the separable approximation. The reason for 
this behavior can be traced to the three-momentum form 
factors commonly used in separable interactions. Stan- 
dard analytic continuation of the KN momentum vari- 
able into the subthreshold region lets these form factors 
increase beyond their threshold (zero-momentum) mag- 
nitudes, thereby enhancing the subthreshold amplitude 
artificially. While such analytic continuations need not 
be performed in the Faddeev approach, this example nev- 
ertheless demonstrates that extrapolations into the far- 
subthreshold region are confronted with off-shell uncer- 
tainties which severely restrict the predictive power of 
KNN binding energy calculations. 

Next, we examine the dispersive effect induced by the 
imaginary part of the KN potential. The advantage of 
the Faddeev method is its ability to deal consistently to 
all orders with the imaginary parts of the interactions 
whereas the variational method can handle this only per- 
turbatively. In Section III C we have estimated that dis- 
persive corrections from these imaginary parts amount 
to about 6 ± 3 MeV binding per nucleon, based on the 
analysis of the two-body KN channel. Thus the disper- 
sive correction to the ppK~ system would add another 
AB(ppK~) < 15 MeV to the total binding energy. 

On top of these effects, there are several minor differ- 
ences between two schemes such as the NN potential, 
admixture of the 1$-) component, and so on. These fac- 
tors might work together to contribute to the difference 
between the present result and that in Ref. [20]. 

We add a short comment on the sensitivity of the ppK~ 
results to details of the NN potentials. In fact, as long 
as the KNN system is only weakly bound, the depen- 
dence on different types of NN interactions is weak. We 
have performed test calculations using a soft-core poten- 
tial quite different from Avl8 but equivalent with respect 
to reproducing low-energy NN data. The resulting prop- 
erties of the ppK" quasibound system turn out to be very 
similar to those obtained with the Avl8 interaction. De- 
tails are given in Appendix D. 



1 This is equivalent to a dispersion relation approach, absorbing 
uncontrolled high-energy behavior in a few subtraction constants. 
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FIG. 3: (Color online) (a) NN density in ppK~ , as function of NN relative distance, calculated with the "Type I" ansatz. 
"ORB", "HNJH" and "BNW" results are shown as crossed, solid and dashed lines, respectively. For comparison, the deuteron 
density calculated with Avl8 potential is depicted as dashed line with filled circle, (b) display of r 2 p(r) for the same densities. 
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FIG. 4: (Color online) KN density in ppK~ , as function of KN relative distance, calculated with the "Type I" ansatz. (a) 
Separate display of each isospin component (I — and 1 = 1). Lines without (with) symbols show the I = (/ = 1) component 
of the KN density. Models "ORB" , "HNJH" and "BNW" are depicted with crossed, solid and dashed lines, respectively, (b) 
Normalized 1 = KN relative density in ppK~ for the model "HNJH" with the "Type I" ansatz (solid line). Solid line with 
diamond shows KN density in A(1405) calculated with the same model. 



D. Corrections from p-wave KN interactions 

A rough leading-order estimate of the effect of p-wave 
KN interactions can be performed by computing expec- 
tation values of suitably parametrized p-wave K~p and 
K~n potentials. The p-wave KN effective potentials 
used here, 

47>,^) = Vexp[-r 2 /^] V , (36) 

with r = r^ — riv and the reduced KN energy u>, involve 
the energy-dependent p-wave scattering "volumes" C p for 
-ftT^-proton and C n for if^-neutron, where C„ ~ 2C p . 
Their prominent feature is the £(1385) resonance. We 
use a parametrization updated from Ref. [41] which has 
also been used in our previous study [22]. The energy 



dependence of C p is shown in Fig. 5. 

Expectation values AV^ N = (*|wj™ vc + v^™ e \^) 

of this p-wave KN potential are then computed with 
the ppK~ wave functions of all four chiral models under 
consideration, changing the range parameter a p within a 
reasonable interval, 0.4 fm to 0.9 fm. Then Re AVgjy is 
found to be small and repulsive, 

1.5 < ReAVgN £ 5.0 MeV. (37) 

The weak binding of the ppK~ system places the effec- 
tive yfs for the KN subsystem above the position of the 
£(1385), hence the positive sign of Re AVg N . The con- 
tribution of p-wave KN interactions is small if the ppK~ 
binding from the driving s-wave interactions is relatively 
weak as it turns out in the present work. If the binding 
were strong enough to move the effective energy of the 
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FIG. 5: Energy dependence of if ~-proton scattering volume, 
C p . Its real (imaginary) part is drawn with a solid line (dotted 
line with diamond). 



KN pair down below the £(1385) resonance, the p-wave 
potential would act attractively and would tend to in- 
crease the binding. The present results appear to rule 
out this possibility. 

For the correction to the ppK~ width from p-wave KN 
interactions one finds Ar = -2 lmAV RN <~ 10 - 35 MeV. 
This relatively large Ar (especially for the type I case) 
results from the fact that the two-body energy of the KN 
subsystem is located close to the £(1385) resonance with 
its prominent imaginary part, see Fig. 5. Obviously we 
can give here only a rough estimate for orientation. 



E. Estimate of antikaon absorption by the 
two-nucleon pair 



A further point of interest is the contribution to the 
decay width of the ppK~ cluster from the two-body ab- 
sorption process K~pp — > YN. This effect is not in- 
cluded in the imaginary part of the KN potential used 
to estimate the decay width (22) in previous sections. 

For first guidance, let us start with the formula of the 
decay width for K~ absorption on proton pairs in a heavy 
nucleus [42-44] 



Ar„»,= = 



2nB 



P pp (lo) j d 3 r PR {r) p 2 N (r), (38) 



where p^ and p R are the one-body densities of nucleon 
and antikaon, u> = ra R — B is the energy of the meson, 
and (3 pp (uj) is a kinematical factor normalized to unity 
at threshold, u = m R [42, 44]. This factor reflects the 
phase space and kinematics for the relevant decay chan- 



nels (£ + n, £°p, and Ap in the present case): 

r (,a - v aY ( M( - mK y 

y=A,s v v 1 

^/[M 2 (uj) -m A N + M 2 ] 2 - 4M 2 (u)M^ 
X y/[M 2 (m K ) -m 2 N + M^] 2 - 4M 2 (m K )M£ 

where M(u) = 2M N + uj and a A = 1, a E = 2. Eq. (38) 
derives from an effective contact interaction, with the 
coupling strength Bo estimated empirically as Bo ~ 1 
fm 4 [44] on the basis of kaonic atom data summarized 
and discussed in Refs. [42, 43]. 

The absorptive width in the form (38) is proportional 
to the probability of finding three particles (antikaon and 
two nucleons) at the same space point, as expressed by 
the product of the one-body densities. This treatment is 
justified for large nuclei where the independent particle 
picture works reasonably well. However, for the few-body 
K~pp system, the following modifications are required: 

(i) correlations between the two nucleons must be 
taken into account, and 

(ii) finite range effects in the absorption process must 
be considered. 

The first modification is mandatory because of the re- 
pulsive core in the nucleon-nucleon interaction. Once 
the correlations are taken into account, a local absorp- 
tive contact interaction (38) is not appropriate because 
it provides almost no width: the probability for two nu- 
cleons to be found at the same spot is basically zero. A 
proper treatment must therefore deal with finite range 
effects in the absorption process which in turn requires a 
more detailed assessment of the underlying microscopic 
mechanisms. 

By analogy with pion absorption in the deuteron, a 
leading microscopic process would be one-pion exchange 
(Fig. 6, left). This process, however, does not contribute 
to K~ absorption on a proton-proton pair with spin 
S = in the K~pp system considered here. So the 
driving absorption mechanism is expected to come from 
exchanges of two pseudoscalar mesons as illustrated in 
Fig. 6. In chiral effective field theory, these are subleading 
one-loop terms in the two-baryon system involving ^7T7t 
or KKtt couplings to one of the baryons and exchange 
of the (interacting) two mesons with the second baryon. 
These processes occur when the external antikaon over- 
laps with one of the nucleons. The effective range of the 
absorption process is then related to the mass spectrum 
representing the exchanged two-meson system. 

Thus, in order to estimate the width from the two- 
nucleon absorption process, we need to extend the for- 
mula (38) to satisfy the requirements (i) and (ii). The 
result should reduce to the original form (38) in the limit 
of no NN correlations and zero-range interaction, so that 
contacts to phenomenology can be made concerning the 
absorptive coupling strength Bo- 
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FIG. 6: Feynman diagrams for typical microscopic mecha- 
nisms of the two-nucleon absorption process, (a): one pion 
exchange which is forbidden in the present system, (b) and 
(c): two-meson exchange. 



First note that the product of one-body densities can 
be expressed as 

Pi?( r ) Pjv( r ) = J d 3 x x J d 3 x 2 p R (r) p N {x x ) p N (x 2 ) 

x 6 3 ( Xl - x 2 ) £ 53{Xi ~ r) . (39) 

i=l,2 



For the modification requested by item (i), we introduce 
the three-body density as 

p ( - 3) (r,x 1 ,x 2 ) 
= < * \S 3 (r R - r)6 3 ( ri - Xl )5 3 (r 2 - x 2 )\ *) 
=p R (r)p N (x 1 )p N (x 2 )[l - C(r,xi,x 2 )] , 

where the function C(r,ri,r 2 ) represents the correla- 
tions among the particles. The product of one-body den- 
sities in Eq. (39) can now be replaced by the three-body 
density p^ (r, Xi,x 2 ) which reduces to the original form 
in the limit of C — > 0. The modification (ii) is imple- 
mented by replacing the delta function in Eq. (39) by a 
finite range distribution, such as a normalized Gaussian: 



d 3 (x 1 - x 2 ) -> G(x 1 - x 2 ; a) 



1 



2 / 2 

/a 



Taking a — ► 0, this distribution turns into the delta func- 
tion. In summary, the formula for the absorptive width 
of the ppK~ few-body system is given by 

AT abs (K-pp -» YN) = ^^/3 pp (u;) 

x J d 3 r J d 3 x p {3) {r,r,x)G{x-r;a) . (40) 

The variation of Ar a b s with respect to changes of the 
Gaussian range parameter a is shown in Fig. 7. Note 
that, from the point of view of the microscopic hit or irK 
exchange mechanisms discussed previously, the charac- 
teristic ranges are covering a band from a ~ 0.4 fm for 
K* exchange to about s ~ 1.4 fm for uncorrelated 2ir 
exchange. For a typical choice of the absorptive coupling 
strength, Bq = 1 fm and a = 0.6 fm, one finds: 




0.8 1.2 
a [fin] 



FIG. 7: (Color online) Absorptive width Ar a {, s for ppK — > 
YN as function of Gaussian range parameter a. Results are 
shown for the models "ORB", "HNJH", "BNW" and "BMN" 
discussed in the text. 



for all models (with Type I or Type II ansatz). The 
coupling strength Bq is, of course, subject to consider- 
able uncertainties. From the recent survey [43] based on 
kaonic atom data analysis one estimates a possible band 
width 0.85 < Bq < 1.5 fm 4 [44] which translates into 



4MeV< Ar afcs < 12MeV 



(42) 



5 MeV < AT abs < 8 MeV 



(41) 



Given this estimated upper limit for Ar a f, s in the ppK~ 
system we conclude that the absorptive width in such a 
dilute, weakly bound system is expected to be less impor- 
tant than the one from the mesonic KNN — > irYN decay 
under the given conditions. This situation may change 
qualitatively for heavier nuclei with potentially stronger 
K binding where the absorptive width can be far more 
prominent as pointed out in Ref. [44]. 



V. SUMMARY 

We have investigated the prototype of antikaonic nu- 
clei, ppK~, using a variational approach with an effec- 
tive KN interaction derived from chiral SU(3) coupled- 
channel dynamics. Several versions of such KN inter- 
actions have been employed in the calculations. They 
all satisfy the necessary constraint of reproducing empir- 
ical K~p threshold information and the 7r£ mass spec- 
trum in the region of the A(1405), within (admittedly 
large) uncertainties of the existing experimental data 
base. Furthermore, a realistic nucleon-nucleon interac- 
tion (Argonne vl8) has been used throughout. This in- 
teraction properly accounts for repulsive short-distance 
NN correlations. 

The primary problem faced in the theoretical part of 
the quest for antikaon-nuclear quasibound systems is the 
subthreshold extrapolation of the s-wave KN interac- 
tion. As a consequence of the strong KN «-> 7r£ coupled- 
channel dynamics, this interaction is complex, non-local 
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and energy dependent. While the driving interaction ker- 
nel is determined by chiral SU(3) dynamics, the off-shell 
extrapolation into the far-subthreshold region is subject 
to uncertainties which limit the predictive power of the 
theory. A minimal condition for any such calculation 
is to account in detail for the coupled-channel dynamics 
that governs the formation of the A(1405) as a quasi- 
bound KN state embedded in the strongly interacting, 
resonant 7rS continuum. When this is done, the resulting 
KN effective potential turns out to be significantly less 
attractive than anticipated in a simple phenomenological 
approach using a local, energy-independent potential. As 
a consequence we arrive at weaker ppK~ binding than 
that previously suggested. 

As an independent test of the variational method ap- 
plied in the present investigation, we have performed cal- 
culations using the phenomenological Akaishi-Yamazaki 
(AY) potential. The stronger binding found with this en- 
ergy independent local potential is indeed reproduced as 
reported in detail in Appendix C. 

The results of the present variational calculation are 
summarized as follows: 

• The calculated binding energy of the ppK~ cluster, 
based on the leading s-wave KN interaction only, 
is 

B{ppK~) ~ 20 ±3 MeV , 

where the error indicates variations using four dif- 
ferent versions of chiral SU(3) coupled-channel cal- 
culations. The decay width into irYN final states 
is estimated to be in the range 

T(ppK- -» ttYN) - 40 - 70 MeV . 



• Differences between these variational results and 
those obtained from Faddeev calculations [20] pre- 
sumably relate in large part to irT,N three-body dy- 
namics not incorporated in the present framework. 
Some part of these differences may be attributed 
to the use of separable approximations for the KN 
and nY interactions when solving Faddeev equa- 
tions, and to dispersive effects not covered by the 
variational approach. We estimate the dispersive 
correction to the total ppK~ binding energy to be 
AB(ppK-) < 15 MeV. 

• This additional binding is partly compensated by 
corrections from p-wave KN interactions which 
have a repulsive effect as long as the effective two- 
body energy in the KN subsystem stays above the 
£(1385) resonance. The overall binding energy, 
after corrections and with conservative error as- 
signment, is then expected to be in the interval 
B{ppK-) - 20 - 40 MeV. 

• The p-wave interactions tend to increase the 
mesonic decay width for ppK~ — ► irYN by an 
amount of 10 - 35 MeV. 



• The additional effects of KNN — > YN absorption 
are estimated to increase the decay width of the 
quasibound ppK~ state further by an mount of or- 
der 10 MeV. 

• In the variational approach the wave function of the 
ppK~ quasibound state can be computed and ana- 
lyzed. Given its weak binding, the system is rather 
dilute. The average distance between the two nu- 
cleons is about 2.1 — 2.2 fm. The short-range re- 
pulsion in the NN interaction plays an important 
role in keeping the nucleons apart. The antikaon 
likes to minimize its distance from either nucleon. 
The isospin 1 = KN density distribution within 
the KNN cluster is reminiscent of the A(1405) in 
the KN two-body system, indicating a A(1405)A 
dibaryonic hybrid structure. However, with its es- 
timated lifetime of order r = 1/r ~ 2 frn/c, this 
structure can exist only over a very short time in- 
terval. 

Taking all theoretical uncertainties into account, we ar- 
rive at the conclusion that narrow, deeply bound ppK~ 
clusters are unlikely to exist. Our calculation predicts 
such systems to be weakly bound and short-lived. With 
a rather low binding energy and a total width between 
60 and 120 MeV, such structures would indeed be diffi- 
cult to detect. One should keep in mind, however, that 
the present calculations rely entirely on the constraints 
provided by the presently available sets of K~p threshold 
data and the poorly known 7r£ mass spectrum. Improve- 
ments of this data base will sharpen the theoretical input 
conditions for deriving the effective subthreshold KN in- 
teraction. One can look forward to further developments 
along these lines as the experimental searches proceed 
with higher precision. 
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APPENDIX A: NN POTENTIALS 

For practical convenience in performing the variational 
calculations, the original Avl8 potential has been repre- 
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FIG. 8: (Color online) Central plus spin-spin NN potentials 
in singlet-even (a) and singlet-odd (b) channels. Solid line: 
original Argonne vl8 potential ("Original"). Points represent 
fits to the original potential with four Gaussians ( "Fitted" ) . 



TABLE V: Best fit parameters for the expansion of the Avf 8 
potential into Gaussians. Strengths are given in MeV, ranges 
in fm. 
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APPENDIX B: CHIRAL COUPLED-CHANNEL 
APPROACH AND REGULARIZATION 
DEPENDENCE 

Here we briefly summarize basics of the chiral SU(3) 
coupled-channel approach [26]. In particular, we study 
the regularization dependence of the KN and 7tE am- 
plitudes, an issue that has repeatedly been raised in the 
literature. 

The starting point is the coupled-channel mcson- 
baryon scattering equation 



FIG. 9: (Color online) L 2 potential in NN potential. Figs, 
(a) and (b) show 1 E and 1 channels, respectively. Solid line 
is the original Argonne vl8 potential ("Original"). Points 
represent the potential used in the present study, which fit 
the original potential with three Gaussians ( "Fitted" ) . 

sented by a series of Gaussians. 

The spin state of the two nucleons is restricted to S = 
in our model, so the expectation value of the operator 
&i ■ &2 in the spin-spin potential is equal to —3. We have 
fitted the radial parts of the central potential plus the 
spin-spin potential, v c x (r) — 3i>jp (r) , with four Gaussians 
of different ranges; the radial part of L 2 potential, v^ 2 (r) , 
is fitted with three Gaussians: 

^(r)-34 5 M = ^< + ^e-^< +SS ' X ) 2 ,(Al) 

71=1 

v L x 2 (r) = £^e-^ 2 ' X > 2 , (A2) 

n=l 

with strengths v <^+ ss > x and v^ 2 ' X and ranges r^ SS ' X 
and r^ 2 ' X , where X indicates the channels 1 E, 1 0. The 
values of the parameters are listed in Table V. For illus- 
tration of the fit quality we draw the potentials in Figs. 8 
and 9. 



Tij — Vij + Vik Gk Tkj, (Bl) 

with the interaction kernel Vij derived from the chiral 
SU(3) meson-baryon Lagrangian, and the loop function 
Gk to be discussed in detail later. In order to derive 
the single-channel KN potential, an effective interaction 
kernel V cS is constructed such that the full solution of 
Eq. (Bl) in the KN channel is reproduced by solving a 
single-channel scattering equation with V cS : 

Tu = V eS + V eS G 1 T n , (B2) 

where % = j = 1 now represents the KN channel. This 
V eS is used in Eq. (26) to determine the strength of the 
KN potential. 

In Rcf. [26], dimensional regularization was used in 
the calculation of the loop function G. Dimensional reg- 
ularization has the advantage that the analyticity of the 
loop function is compatible with the dispersion relation 
used in the N/D method [29]. It is however instructive to 
study the regularization dependence in this framework, in 
order to understand possible differences between present 
results for the K~pp system and the results of the Fad- 
deev calculation with interactions constrained by chiral 
dynamics [20] . An analysis of the regularization depen- 
dence in this framework can be found in Ref. [45] . Note 
that differences in the regularization schemes only change 
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the loop function G and leave the interaction kernel V 
untouched. 

Using dimensional regularization, the loop function G 
is given by 

G = — la M + In — + - In — 



g ln ^++(s)</ ) +-(g) 1 
\/s (l>-+(s)<f> — (a) [' 



(B3) 



where a(/i) are subtraction constants, /i is the renormal- 
ization scale, and we have defined 

A = M 2 — m 2 , 



(s) = ±s ± A + 2g\/s , 



*/(s - (M - to) 2 )(s - (M + m) 2 ) 

9 = vi • 

With a sharp three-momentum cutoff, the loop function 
is 

^ <3d . „ 2M f mM A, Af(l + f") 
-ln[(l+r)(l + e M )] 



0™(s)0^(s) J ' 
where g maa; is the three-momentum cutoff and 
«^( S )=± S TA + 2gVir , 



(B4) 



,1/ 



M 2 



A smooth cutoff can be introduced using the Pauli-Villars 
method which corresponds to multiplying a monopole 
form factor 



q 2 - A 2 

to the loop function, leading to 

(47r) 2 [Asm s A 



/i 0_ + (s)0 — (s) 

gA ln </>A,++(s) ^A, + -(g) 
\/s <^A,-+(s) </>A, (s) 



(B5) 



where 



A \ = Af 2 - A 2 , 
0A,±±(«) =±s±A A + 2g A \A 



q\ = 



y/(s - (M - A) 2 )(s - (Af + A) 2 ) 
2^ 
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FIG. 10: (Color online) Comparison of the loop functions: di- 
mensional regularization (solid lines), thee-momentum sharp 
cutoff (dashed lines) , and smooth monopole form factor (dot- 
ted lines), (a): KN channel, (b):, 7rE channel. 



The real parts of the loop functions (B3), (B4), and 
(B5) for KN and 7rE channels are plotted in Fig. 10. 
The imaginary parts are determined by the phase space 
of intermediate meson-baryon states and therefore inde- 
pendent of the regularization procedure. The parameters 
for the dimensional regularization and three momentum 
cutoff schemes are taken from phenomenologically suc- 
cessful models [28, 36]: 



a KN 
?max 



— — 2 , OttS 

= 630 MeV. 



-1. 



fi = 630 MeV , 



The loop functions G dtm and G 3d are evidently quite 
similar. Furthermore, choosing the parameter of the 
monopole form factor as 



A^ =750 MeV, 



A nS = 500 MeV 



the corresponding loop functions G pv behave quan- 
titatively similar as those with dimensional/three- 
momentum cutoff schemes, as seen in Fig. 10. The cutoff 
scales of the form factors (several hundreds of MeV) is 
typical and naturally expected from meson-baryon phe- 
nomenology. 

In summary, we find that smooth cutoff schemes pro- 
vide similar loop functions as those with sharp cutoff or 
dimensional regularization. Since differences in the regu- 
larization schemes only affect the loop functions G in the 
present framework, it is clear that the potentials derived 
in Ref. [26] remain unchanged when adopting a smooth 
cutoff. 

In this respect, the difference between the present re- 
sults for K~pp binding and those in Faddeev method [20] 
cannot be attributed to the regularization method. It is 
worth noting that early studies using the coupled-channel 
approach [27] gave a stronger K subthreshold amplitude 
than the one derived in Ref. [26]. The calculations of 
Ref. [27] used a separable approximation for the potential 
together with the non-relativistic Lippmann-Schwinger 
equation, a feature also shared by the Faddeev calcu- 
lations of Ref. [20]. Thus the larger binding energy of 
K~pp found in Ref. [20] is likely to be related to the type 
of scattering equation and the separable approximation 
used in that work. 
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TABLE VI: Test of the present variational calculation with 
the AY potential [21]. "Model I" and "Model II" are 
the results using the variational wave function described in 
Sec. II A. "ATMS" refers to the results quoted in Ref. [21]. 
B(ppK~) and F are the total ppK~ binding energy and the 
KNN — » nYN decay width. The antikaon binding energy 
Bk is defined in Eq. (16). All the remaining quantities are 
specified as in Table III. All energies are given in units of 
MeV. The r.m.s. distances -Rjvjv and R KN are m f m - 



c 


Model I 



Model II 
finite 


ATMS 


B(ppK~) 


39.0 


51.4 


48 


r 


60.0 


61.0 


61 


B K 


65.8 


80.0 


68 


Tnuc 


46.7 


47.8 






147.0 


162.4 


167 


V{NN) 


-19.8 


-19.2 


-19 


V{KN) 


-166.2 


-194.6 


-196 


Rnn 


1.75 


1.83 


1.90 


Rrn 


1.54 


1.55 


1.57 


P(T N = 0) 





5.9 % 





APPENDIX C: BENCHMARK TEST WITH THE 
AY POTENTIAL 

As a test for the variational method and trial wave 
functions applied in the present work, we perform a cal- 
culation with the KN interaction [the Akaishi-Yamazaki 
(AY) potential] and the NN interaction (the Tama- 
gaki potential) used in Ref. [21]. The AY potential is 
an energy- independent, local potential based on phe- 
nomenology but not constrained by chiral SU(3) dynam- 
ics. 

The results are summarized in Table VI. The first 
two columns ("Model I" and "Model II") show our re- 
sults; the last column ("ATMS") is the original result 
reported in Ref. [21]. In "Model I" the energy variation 
is performed restricting the wave function to its dominant 
|$ + ) component, i.e. with no admixture of the |$_) state 
(the one with NN coupled to isospin T/v = 0), keeping 
the coefficient C = in Eq.(l). "Model II" includes the 
|$_) component, with C determined variationally. In all 
calculations, the convergence of the Gaussian expansions 
(6) and (7) has been checked and found satisfactory with 
N N =N K = 9. 

One evidently finds a high degree of consistency be- 
tween "Model II" and "ATMS" , confirming that the dif- 
ferent variational methods used here and in Ref. [21] are 
of comparable quality. 

The importance of the TV = component of the wave 
function is underlined by the comparison between "Model 
I" with "Model II". Although this admixture is only 




r [fm] 

FIG. 11: (Color online) Comparison of Minnesota potential 
with Avl8 potential. 1 E channel of both potentials are shown. 
Minnesota potential is depicted with a solid line with dia- 
mond. Avl8 potential is depicted with a solid line. 

about 6 % as shown in the last row [P(Tn — 0)] in Table 
VI, switching it off reduces the binding energy by more 
than 20 %. The matrix elements of the KN interaction 
taken between normalized states 1$+) and |$_) are as 
follows: 

($+| v RN |$+) = + i^, (ci) 

(*_| V RN |$_) = \v R ^ + \v R ^, (C2) 

Vkn I*-) - ^ (4/v - 4*) ■ (C3) 

The mixture implied by Eq. (C3) increases the ppK~ 
binding energy. For the AY potential the values of these 
three matrix elements are —173.3 MeV, —112.3 MeV and 
—53.0 MeV, in this order. Taking the mixing ratio be- 
tween and |$_) into account, the actual contribu- 
tions from each matrix element to the KN potential en- 
ergy are -163.1 MeV, -6.6 MeV and -24.9 MeV, re- 
spectively, and so the coupling matrix element (C3) is 
found to be attractive and non-negligible. 

In concluding this Appendix we note again that the 
overall attraction produced by the V RN based on chiral 
SU(3) dynamics, and used in the present work, is consid- 
erably weaker in comparison and leads to &ppK~ binding 
energy less than half of that found with the simple AY 
potential. 

APPENDIX D: DEPENDENCE ON THE NN 
POTENTIAL 

Test calculations have been performed replacing the 
Avl8 by a soft-core NN potential, the Minnesota po- 
tential [46] (see Fig. 11). This potential reproduces low- 
energy JVJV data (scattering lengths, effective ranges, and 
deuteron properties) just like the more realistic Avl8 in- 
teraction. In the actual computations, we have fixed a 



17 



TABLE VII: Dependence of ppK~ results on the NN poten- 
tial. The fist and second lines indicate the type of KN and 
NN potential, respectively. "Minnesota" in the second line 
refers to the soft-core potential of Ref. [46]. E nuc is the to- 
tal energy of nuclear part, namely E nuc — T nuc + V(NN). 
Other quantities shown in this table are the same as those in 
Table VI. All results are obtained with Gaussian expansions 



N N = N K = 


= 5 in Eqs. 


(6) and (7). 
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parameter of the Minnesota potential (u = 1) so that it 
corresponds to a Serber-type potential. 



Results of ppK calculations are summarized in Ta- 
ble VII for the chiral KN potential (HNJH) and for the 
phenomenological AY potential. The sensitivity to de- 
tails of the NN interaction turns out to be marginal for 
the weakly bound ppK~ system and slightly more pro- 
nounced but still weak for the more strongly bound AY 
case. The qualitative difference between strong and soft 
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r [fm] 

FIG. 12: (Color online) NN density for Avl8 potential (solid 
line) and Minnesota potential (solid line with diamond). A 
chiral-based KN potential (HNJH) is used. 

short-range repulsive core becomes apparent, however, 
when examining the NN density distribution within the 
KNN clusters (see Fig. 12). 
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